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Quadratic algebras of the type [Q , Q±] = ±Q±, [Q+, Q-] = aQl + bQo + c 
are studied using three-mode bosonic realizations. Matrix representations 
and single variable differential operator realizations are obtained. Examples 
of physical relevance of such algebras are given. 



1 Introduction 

In recent times there has been a great deal of interest in non-linear defor- 
mations of Lie algebras because of their significant applications in several 
branches of physics. This is largely based on the realization that the phys- 
ical operators relevant for defining the dynamical algebra of a system need 
not be closed under a linear (Lie) algebra, but might obey a nonlinear alge- 
bra. Such nonlinear algebras are, in general, characterized by commutation 
relations of the form 

[T l) T J }=C l] {T k ] , (1.1) 

where the functions CV,- of the generators {T k } are constrained by the Jacobi 
identity 

[T h C jk ] + ft, C ki ] + [T k , Cy = . (1.2) 

The functions can be infinite power series in {T^} as is in the case of 
quantum algebras (with further Hopf algebraic restrictions) and g-oscillator 
algebras. When {C^} are polynomials of the generators one gets the socalled 



polynomially nonlinear, or simply polynomial, algebras. A special case of 
interest is when the commutation relations ( |1 . 1| ) take the form 

[T if Tj] = c%T k , [Tj, T a ] = tlT p , [T a , Tp] = f a p (T k ) , (1.3) 

containing a linear subalgebra. Simplest examples of such algebras occur 
when one gets 

[N , N±) = ±N ± , [N+, N_] = f (N ) . (1.4) 
In general, the Casimir operator of this algebra (|1.4j) is seen to be given by 
C = N + N. + g(N -l)=N-N+ + g(N ) , (1.5) 



where g (No) can be determined from the relation 



g(No)-g(No-l) = f(N ) . (1.6) 

If / (N ) is quadratic in N we have a quadratic algebra and if / (N ) is 
cubic in N we have a cubic algebra. The nonlinear algebras, in particular 
the quadratic and cubic algebras, and their representations have been stud- 
ied in connection with several problems in quantum mechanics, statistical 
physics, field theory, Yang-Mills type gauge theories, two-dimensional inte- 
grable systems, etc. (JTJ-0)- ^ s * n ^ ne case °^ ^ ne ^ e a lg e bras, compact 



nonlinear algebras have finite dimensional unitary irreducible representations 
and noncompact algebras have infinite dimensional unitary irreducible repre- 
sentations. These nonlinear algebras have also been found to be very useful 
in quantum optics with the observation that quantum optical Hamiltonians 
describing multiphoton processes have symmetries described by polynomially 
deformed su(2) and su(l, 1) algebras [H]. Coherent states of different kinds 



of nonlinear oscillator algebras have been presented by several authors (|J 
[p0|1). Recently, a general unified approach for finding the coherent states of 
polynomial algebras, relevant for various multiphoton processes in quantum 



optics, has been presented by us [24 



Algebras of the type ( |1.4| ) have been studied in mathematics literature 



also and shown to have a surprisingly rich theory of representations [p5 
Algebras of the type ( |1.4|) with commutator [N + ,N_] replaced by the anti- 
commutator {iV + , N_} leading to polynomial deformations of the superalge- 
bra osp(l\2) have also been investigated ||26j| . As is well known, in the case 
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of classical Lie algebras bosonic realizations play a very useful role in the 
representation theory and applications to physical problems. The main pur- 
pose of this work is to study some aspects of quadratic algebras relating to 
three-mode bosonic realizations, and the associated matrix representations, 
differential operator realizations, coherent states, and physical applications. 
In Section 1 we briefly review the two-mode bosonic construction of su{2) and 
su(l, 1) algebras. Following closely our earlier work p4| , in Sections 2 and 3 
we construct the compact and noncompact quadratic algebras, respectively, 
using three bosonic modes and study the associated representations. In Sec- 
tion 4 we briefly discuss the coherent states of these algebras. Finally, in 
Section 5 we conclude with a few remarks on the physical and mathematical 
significance of these quadratic and higher order nonlinear algebras. 



2 Two-mode construction of su(2) and 
su(lj 1) : A brief review 

Let us briefly recall the study of su(2) and su(l, 1) in terms of two-mode 
bosonic realizations, to fix the framework and notations for our work. The 
su(2) algebra, 

[Jo,J±]=±J±, [J + ,J-] = 2J , (2.1) 
is satisfied by the Jordan- Schwinger realization 

Jo — - (a\ai — a^a^j , J+ = a[a 2 , J- = J+ = Oia| . (2.2) 
In this realization the Casimir operator becomes 

C = J 2 = J + J- + Jo (Jo — 1) = ^ {a\a\ + a^a^ {a\a\ + a 2 a 2 + 2^ . (2.3) 

Consequently, the application of the realization ( |2.2| ) on a set of 2j + 1 two- 
mode Fock states |n 1; n 2 ) with constant n\ + n 2 = 2j leads to the (2j + 
l)-dimensional unitary irreducible representation for each j = 0, 1/2, 1, . . . . 
Thus, with {\j, m) = \ni = j + m, n 2 = j — m) \ m = j, j — 1, . . . , — j} as the 
basis states, one gets the j-th unitary irreducible representation 

Jo I j, m) = m | j, m) , 

J± I j, m) = y/(jTm)(j ±m + 1) \j, mil), 

™> = j,j-l,... , -j , (2.4) 
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associated with the constant value j(j + 1) for the Casimir operator J 2 . 

Let us now consider the single variable differential operator realization 
corresponding to the above matrix representation ( |2.4p . With the Fock- 
Bargmann correspondence 

(at.a) — ^), |n> — (2.5) 

we can make the association 

J+m j~m 2j i I \j+m 

Z x Z 2 Z 2 \Z\jZi) 



\j, m) 



<J(j + m)\(j - m)\ yJ(J + m)\(j - m)\ 

m = -j, + ,j- (2.6) 



Since j is a constant for a given representation we can rewrite the above as 
a mapping to monomials 

z n 

|j,m)— ^, w (z)= n = 0,1,2,... ,2j. (2.7) 

v /n!(2j-n)! 

Then, it is obvious that the above set of (2j + 1) monomials ( |2.7| ) forms the 
basis carrying the finite dimensional representation ( |2.4| ) corresponding to 
the single variable realization 

J = z—-j, J + = -z 2 ^- + 2jz, J- = -f- (2.8) 

(JL/C KAJ/C Ul/C 

In an analogous way, for su(l, 1) the two-mode bosonic realization 
K = — (a\ai + a\a 2 + l) , K + = a\a 2 , K_ = K\_ = a\a 2 , (2.9) 



2 

satisfies the algebra 

[K , K±) = ±K± , [K+, K.] = -2K . (2.10) 



The Casimir operator is 

C = K 2 = K + K^ - K (K - 1) = i 



(2.11) 
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Consequently, the application of the realization ( |2.9| ) on any infinite 
set of two-mode Fock states {| k, n) = \n, n + 2k — 1)| n = 0, 1, 2, . . .} (or 
{| k, n) = \n + 2k — 1, n)\ n = 0, 1, 2, . . .}), with constant \n\ — n 2 \ = 2k — 1, 
leads to the infinite dimensional unitary irreducible representation, the so- 
called positive discrete representation T> + {k) , for each k = 1/2, 1, 3/2, . . . . : 

K \k,n) = (k + n)\k,n) , 



K+\k,n) = yj(2k + n)(n + l) \k,n + l) , 

K-\k,n) = \]{2k + n- l)n \k,n- 1) , 

n = 0,1,2,.... (2.12) 

with the Casimir operator K 2 taking the value k(l — k) in the representation. 
As in the su(2) case, we can make the association 

r n r n+2fc-l / \n 2fc-l 

|ib,n)— > Z 1 ^ 2 = ; ( "^ 2) " 2 (2.13) 

fn\(n + 2k-l)\ Jn\(n + 2k-l)\ 



Then, with k being constant in a given representation, it is obvious that the 
infinite set of monomials 

ijj Kn {z) = -= n = 0,l,2, (2.14) 

yfn\(n + 2k-l)\ 



forms the basis carrying the representation (|2.12|) corresponding to the single 
variable realization 

K = z^- + k, K + = z, K. = z^— + 2k4- ■ (2.15) 

az az z az 

3 Three-mode construction of compact 
quadratic algebras 

A quadratic algebra is defined, in general, by the commutation relations 

[Qo, Q±] = ±Q± , [Q+, Q-] = aQ 2 + bQ + c, (3.1) 
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where (a, b, c) are constants, or commute among themselves and with 
(Qo, Q+, Q-) so that they take constant values in any irreducible representa- 
tion. When a = c = and b = ±2 one gets, respectively, su{2) and su(l, 1) 
as special cases. Let us now construct a class of compact quadratic algebras 
using three bosonic modes. To this end, let 

Qo = - (a\ai + a\a 2 - 2a\a 3 + l) = - (K - 4a 3 J , 
Q + = a\a} 2 a 3 = K + a 3 , Q- = Q+ = aia 2 al = K_a\ , 




C = - (a\a x + a\a 2 + 2a\a 3 + l) = - (k q + a|a 3 ) , (3.2) 

where (K ,K + ,K_) generate su(l, 1) with K 2 as the Casimir operator. 
Then, it is found that 

[£,£] = 0, [£,Q 0> ±] = 0, [£,Q ,±] = 0, 
[Qo,Q±] = ±Q±, [g+,g_]=3Qg + (2£-l)Q +(/C-£(£ + l)). 

(3.3) 

Thus, (Qo, Q+; Q-) generate a quadratic algebra with K and £ taking con- 
stant values in any irreducible representation. The Casimir operator of this 
algebra is given by 

C = Q+Q- +Ql + {£- 2)Q 2 + (JC - C 2 - 2C + l) Q , (3.4) 

apart from an additional constant function of K, and C, following the recipe 
contained in ([L4D-(|OD. 

The condition that K, and L take constant values in an irreducible repre- 
sentation fixes the basis to be the set of three-mode Fock states 

\k,l,n) = \n,n + 2k-l,2l-k-n) , n = 0, 1, 2, . . . , (21 - k) , (3.5) 

with 

2Z-fc = 0,1,2, ... , k = 1/2,1,3/2, ... , (3.6) 

and 

JC \k, I, n) = k(l — k) \k,l,n) , C \k, l,n) = I \k, I, n) . (3.7) 
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The basis states ( |3.5| ) carry the (21 — k + l)-dimensional unitary irreducible 
representation of the quadratic algebra ( |3.3| ) which can be labeled by the 
values of the pair (k,l). Explicitly, the (k,l)-th representation is: 



Qo\k,l,n) = (k — I + n) \k,l,n) , 



Q + \k,l,n) = \J(n + l)(n + 2k) (21 - n - k) \k,l,n + l) , 



Q-\k,l,n) = yjn(n + 2k - l)(2l - n - k + 1) \k, l,n-l) . 
fC\k,l,n) = k(l — k) \k, I, n) , C\k,l,n) = l\k,l,n) , 

n = 0, 1,2, ... ,(21 -k). 



(3.8) 



The Casimir operator has the value (/ 3 + (I + 1) [Ar(l — k) — 1] + 1) in this 
representation. 

It should be noted that for k > 1/2 the same representation is obtained 
by the choice of basis states as \k, I, n) = \n + 2k — 1, n, 21 — k — n) , with 
n = 0, 1, 2, ... , (21 — k) . It should also be noted that, unlike in the case of the 
su(2) algebra, the quadratic algebra (|3.3l ) has infinitely many inequivalent 
unitary irreducible representations of the same dimension. For example, for 
each value of k = 1/2, 1, 3/2, . . . , there is a 2-dimensional representation of 
the algebra (|3.3|) given by 




k(l-k), C = l 




(3.9) 



as can be verified directly. 



As before, let us make the association 



k, I, n) 



2k-\Jll-k 
2 z 3 



(z x z 2 lz z ) 



n 



(3.10) 




Since k and / are constants for a given representation we can take 



z 



ii 





n= 0, 1,2, ... ,(21- k), 



(3.11) 
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as the set of basis functions for the single variable realization 
Qo = z^ + k-l, Q + = -z 2 -^- + (2l-k)z, Q_ = z ^ + 2k^, (3.12) 



leading to the representation 

4 Three-mode construction of noncompact 
quadratic algebras 

Let us choose 

Qo = j (a\ a i + a l a 2 + 2a\a 3 + l) = - (k + 4 a s) , 

Q + = ala^al = K + a 3 , Q_ = Q+ = a\a 2 a 3 = K_a 3 , 

K. — — 1 — (a\ai — a\a<2^ = K 2 , 

£ = - (alat + a\a 2 - 2a\a z + l) = - (K - a\a^ , (4.1) 

where (K Q , K + , KJ) generate su(l, 1) with K 2 as the Casimir operator. Now, 
the algebra becomes 

[K,£] = 0, [K,Q , ± ] = 0, [£,Q 0) ±] = 0, 
[Qo,Q±] = ±Q±, [Q + ,Q-] = -3Ql-(2£+l)Q -(IC-£(£-l)). 

(4.2) 

The Casimir operator of this algebra is given by 

C = Q+Q_ -Q 3 + (£- l)Ql - (JC - £ 2 ) Qo (4.3) 
apart from additional constant functions of K, and £. 

The constancy of /C and £ in an irreducible representation leads to the 
choice of basis consisting of the infinite three-mode Fock states 

\k,l,n) = \n,n + 2k-l,n + k-2l) , n = 0,1,2,..., (4.4) 

with 

k-2l = 0,1,2,... , k = 1/2, 1,3/2,... , (4.5) 
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and 

K \k, I, n) = k(l - k) \k, I, n) , C \k, I, n) = I \k, I, n) . (4.6) 
The set of basis states Q4.4|) carry the infinite dimensional unitary irreducible 



representation of the quadratic algebra ( |4.2|) given by 
Qo\k,l,n) = (k — I + n) \k,l,n) , 



Q + \k,l,n) = ^J(n+l)(n + 2k)(n + k-2l + l) \k,l,n+l) 

Q_\k,l,n) = yjn(n + 2k-l)(n + k- 21) \k, l,n-l) , 
K \k, I, n) = k(l — k) \k, I, n) , C\k,l,n) — l\k,l,n) , 

n = 0,l,2,.... (4.7) 

The Casimir operator has the value / (I — k 2 ) in this (k, Z)-th representa- 
tion. As in the compact case, to obtain this representation one can also use 
an alternative set of basis states \k, I, n) = |n, n + 2k — 1, n + k — 21) , with 
n = 0,1,2, .... 

From the association 

\k,l,n) — > , = (4.8) 



^n\(n + 2k-l)\(n + k-2l)\ 



it is clear that we can take 



1pk,l,n( z ) = 1 = ( 4 - 9 ) 

^Jn\(n + 2k- l)\(n + k-2l)\ 

as the set of basis functions for a single variable realization of the algebra 
( |4.2|) . The corresponding realization is 

Qo = z-^ + k-l, Q + = z, 

dz 

Q- = z 2 ^ + (3k-2l + 2)z-^+(2k 2 -Akl + 2k)-^-. (4.10) 

dz 4 dz 1 v 7 dz 



5 Coherent states of the quadratic algebras 

As is well known, coherent states form an overcomplete set of nonorthogonal 
states, generally labeled by a continuous index, say a, providing a resolution 
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of identity / |a) (a) da = I and hence a useful set of basis states. Coherent 
states associated with the unitary irreducible representations of the dynam- 
ical algebra of a physical system are very useful in certain studies of the 
system. 

For the noncompact algebra (|4.2j ) the Barut-Girardello type coherent 
states associated with the (k, /)-th irreducible representation can be defined 
as the eigenstates of Q- : 



where the complex a labels these states. Solving for \k, I, a) using the repre- 
sentation ( |4.7| ) we get easily 





(5.1) 



n=0 



k, I, a) 




oo 



a 



n 




k,l,n) . (5.2) 



In terms of the single variable realization the coherent state equation ( |5 . 1| ) 
can be written as 




which is the equation for 



tffc,j(af)*) 



F 2 (-; 2k, k - 21 + 1; az) . 



(5.4) 



The resolution of the identity is given by: 




(5.5) 



With a polar decomposition ansatz a = re 10 we get, 




(5.6) 
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where 

Ni, k (r 2 ) =o F 2 (-, 2k, k-2l + l; r 2 ) (5.7) 
The integral (5.5) reduces to the following condition on a(r 2 ): 

1 f°° 2 , f 2W 2Nfn+ll-l 1 TV i 1 A r ( 2fc + n ) T ( k ~ 2/ + 1 + n) 

' d(r ) cr(r ) (r 2 ) ( " +iJ 1 = — r(n + 1) 



2 7 v ' y ' y ' 2vr ' Y{2k)Y{k - 2/ + 1) 

(5.8) 

<r(r 2 ) is found by an inverse Mellin transform to be: 

a(r2) = vrr(2A;)(r(A: -21 + i) G 03^ 2 |^-2Z + i,2 fc ), (5-9) 



where G is the Meijers G function ||21|| . 



Another useful set of coherent states for the quadratic algebras will be 
the analogues of the well known Perelemov-type coherent states. A technique 
for generalizing the Perelemov-type construction of coherent states for the 
nonlinear algebras, using a mapping of the given compact or noncompact 
algebra to su{2) or su(l, 1) respectively, has been presented in detail in our 
earlier work |p3 |. For completeness we give the expressions for these states. 



For the non compact case the Perelomov type states are: 
|/3 = iV 



r<» +B )r(* -a + i +B)lJ|itB>| 



y\ T(n + l)T(2k)T(k~2l + r 
where the normalization coefficient N is given by 

^=( 2 F (2M-2/ + l;(|/3| 2 )))-i (5.11) 
The resolution of the identity in this case reduces to finding a{r 2 ) such that, 

. d(r 2 ) a(r 2 ) (r 2 )^ 1 = -T(n + 1) - ~ 2* + 1) ( } 

o v ; v M ; 7r ; T(2k + n)T(k - 21 + 1 + n) v ; 

and the resultant a(r 2 ) is given by 

a(r 2 ) = -(T(k -21 + l)T(2k)G 1 2 '° 1 (r 2 \ k - 2l+1 ' 2k - 1 ) } (5.13) 
where G is the Meijers G function. 



11 



The corresponding states for the compact case are given by: 

2l-k 



| of, k, I >= N ( a ) 

n=0 



(21 - k)\(2k -1 + n)! 
^n!(2Z-fc-n)!(2fc-l)! 



\l, k.n > 



(5.14) 



For the purposes of calculating the measure for the resolution of identity we 
define 7 = -. The coherent state \j,k,l > becomes 



2l-k 



\ 1 ,k,l>=N 1 k - 2l J2(l) r 



n=Q 



(2l-k)\{k + 2l-l-n) 
V n\{2l- k-n)\{2k- 1)! 



■|Z,Jfe.n>, (5.15) 



with the normalization coefficient N given by: 
J(k + 2l) 



N 



(l7l 



2\2l~k 



$(A;-2Z,l-2Z-fc;(|7| i )))^. 



(5.16) 



where <&(a,b,x) is the Confluent Hypergeometric function ( 1 i 7 i). 

The resolution of the identity for the coherent states I7, k, I) is given by 



rf/i(7,7*; k,l)\r, k,l)(r, k,l\ = lr k , 



where 1^ is the projection operator on the subspace H.2i-k- 

2l-k 

h,k = \l,k,n)(l,k,n\. 

n=0 

16 



(5.17) 



(5.18) 



Again defining 7 = re we have, 



2 ^ T(2l + k-n) 
t Q n\{2l - k - n)\ 

where we have defined 
M(r 2 ;kJ) 



r 2 ) n M(r 2 ;k,l)d(r z 



\l, k,n)(l, k,n\ = (5.19) 



n(2l-k)\ 



a(r 2 ;k,l) 



T(2l + k) ($(fc-2Z,l-2/-A;;r 2 ))' 



By using the integral ||23fl , 



POO 

/ r^Q^^-^dr 
Jo 



r(6)r(c)r(a - b) 

T{a)T{c-b) ' 



(5.20) 



(5.21) 
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we obtain 



M{r 2 ;k,l) 



T{21 -k + 2) 



<&{2l-k + 2;2l + k + l;-r 2 ). 



T(2l + k + V 

This gives us the final expression for the integration measure: 

d/x( 7 , 7*; k, I) = ^l + kVlf ^' 21, l - 2l ~ k] r2 M 2l - k + 2 'i k+2l+l; -r 2 )d{ 

(5.22) 

The resolution of the identity is important because it allows the use of 
the coherent states as a basis in the state space. All three sets of coherent 
states can also be easily shown to be overcomplete |24j . 



6 Concluding remarks : Physical applications 

As already mentioned in the introduction the nonlinear algebras are finding 
several applications in physical problems. Here we shall like to make some 
observations with reference to the bosonic constructions of the quadratic al- 
gebras we have presented. 

Let us consider the Hamiltonian 

H = a\ai + a\a2 + 20303 (6-1) 

which describes, with units h = 1 and u — 1, and apart from the additional 
zero-point level constant, a three-dimensional anisotropic quantum harmonic 
oscillator with the frequency in the third direction twice that in the perpen- 
dicular plane. From the construction of the compact quadratic algebra given 
above we recognize that 

H = 4C-1, (6.2) 

where (£, /C, Qq, Q±) generate the quadratic algebra (|3.3|) . Thus, (/C, Qq, Q±) 
are constants of motion for the system ( |6.1|) and the quadratic algebra (|3.3| ) 
is its dynamical algebra. Since C has the spectrum 

£ = Z = n /4, 71 = 1,2,3,..., (6.3) 

it is clear that the Hamiltonian (|6.1| ) has the spectrum 

H = N, N = 0,1,2,... . (6.4) 
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Each level can be labeled by the eigenvalues of the complete set of commut- 
ing operators (H, /C, Q$). 

It is interesting to compute the degeneracy of the iV-th level using the 
representation theory of the algebra ( |3.3| ). For the iV-th level the value of £ 
is / = (N + l)/4. Calculating the corresponding values of k for which finite 
dimensional representations are possible we find that the dimensions of the 
associated irreducible representations are (1,2,... , 2m + 1) if N — 4m or 
4m + 1 , and (1,2,... , 2m + 2) if N = Am + 2 or Am + 3. The degeneracy of 
the level is the sum of the dimension of the k = 1/2 representation and twice 
the dimensions of k > 1/2 representations. One has to count the dimensions 
of k > 1/2 representations twice in the sum since there are two possible 
choices for the bases leading to the same representation in these cases as 
already noted. Now, the four cases, N = Am, Am + 1,4m + 2,4m + 3, are 
to be considered separately. The result is : the degeneracies of the levels, 
iV = 4m, 4m + 1,4m + 2, 4m + 3, resectively, are (2m + 1) 2 , (2m+ l)(2m + 2), 
4(m + l) 2 and 2(m + l)(2m + 3) . In other words, the number of composi- 
tions of the integer N (partitions with ordering taken into account) in the 
prescribed pattern n\ + i%2 + 2ns, with the interchange of n\ and ri2 taken 
into account, is (2m+l) 2 , (2m + l)(2m + 2), 4(m+l) 2 and 2(m+ l)(2m + 3), 
if N = Am, Am + 1, 4m + 2, and 4m + 3, resectively. 

It is to be noted that in the above example the sum of all the dimensions 
of the irreducible representations associated with the given I = (N + l)/4 
gives the number of partitions of N in the pattern ^1+^2 + 2713, disregarding 
the interchange of ni and n^. This leads to the result that the number of 
such partitions is (m + 1) (2m + 1) for N = 4m or 4m + 1 and (m + 1) (2m + 3) 
for iV = 4m + 2 or 4m + 3. Thus, it is interesting to observe this connection 
between the quadratic algebra and the theory of partitions. 



An another interesting possibility is suggested by the structure of the 
compact quadratic algebra (|3~3|). Let us define 



N = Q , A 

The the algebra 
[N,A] 



1 



£(£ + l)-/C 
becomes 



Q-, ^ 



£(£ + !)-£ 



Q+. (6.5) 



-A 



~N, A^} = A ] 
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A, A* 



1 



C{C + \)-K 



N 



2C - 1 
£{£ + !) -K 



N . 



(6.6) 



We may consider this as the defining algebra of a quadratic oscillator, corre- 
sponding to a special case of the general class of deformed oscillators 

0): 



[N,A] 



-A 



N, A^ =AK A, A* 



F(N) . 



(6.7) 



The quadratic oscillator (|6.6|) belongs to the class of generalized deformed 
parafermions jl2| . It should be interesting to study the physics of assemblies 
of quadratic oscillators. In fact, the canonical fermion, with 



N 




1 



/ 



1 




is a quadratic oscillator! Observe that 

[NJ} = -f, [tf,/tl=/t, \f,P 



1 





1 



1 



-N - -N 



(6.8) 



(6.9) 



To conclude, we find that the study of polynomial algebras should lead 
to interesting connections with the theory of partitions and new physical 
systems in which the polynomial oscillators replace the role of the canonical 
oscillators, besides having the other known applications. 
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